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ABSTRACT 


piis  Note  contains  curves  of  buckling  coefficients  and 
formulas  for  the  calculation  of  the  critical  external 
pressure  of  finite  length,  circular  cylindrical  shells 
with  sandwich  walls.  The  facings  of  the  sandwich  are 
isotropic,  of  equal  or  unequal  thickness,  of  the  same 
or  different  material,  and  their  individual  stiffness 
is  not  taken  into  account.  The  sandwich  core  is  iso¬ 
tropic  or  orthotropic  having  natural  axes  in  the  axial, 
tangential,  and  radial  directions  oi  the  cylinder.  If 
the  cores  are  very  rigid,  the  method  yields  results  that 
are  substantially  those  of  von  Mises. 
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Introduction 


Design  curves  for  the  critical  external  radial  pressure  of  circular  cylin¬ 
drical  shells  with  sandwich  walls,  calculated  according  to  the  formulas 
developed  at  the  Forest  Products  Laboratory  (9)  Jt  are  presented  in  (his 
report.  The  sandwich  cylinder  walls  have  isotropic  facings  of  equal  or 
unequal  thickness,  and  of  the  same  or  different  materials,  and  orthotropic 
or  isotropic  cores.  It  is  assumed  that  Poisson's  ratio  is  the  same  for 
both  facings.  The  natural  axes  of  the  orthotropic  cores  are  axial,  tangen¬ 
tial,  and  radial.  These  formulas  reduce  substantially  to  those  developed 
by  von  Hises  (7,  10,  14)  when  the  core  ~8  very  rigid. 


ijhis  Note  is  the  latest  revision  of  "Design  Curves  for  the  Buckling  of 
Sandwich  Cylinders  of  Finite  Length  under  Uniform  External  Lateral  Pres¬ 
sure,"  by  Charles  B.  Norris  and  John  J.  Zahn.  It  was  originally  issued 
as  Forest  Products  Lab.  Rpt.  1869  in  1959,  and  revised  as  U.S.  Forest 
Service  Research  Note  FPL-07  in  1963. 
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—This  progress  report  is  one  of  a  series  (ANC-23,  Item  57-3)  prepared  and 
distributed  by  the  Forest  Products  Laboratory  under  U.S.  Navy,  Bureau  of 
Naval  Weapons  Order  No.  19-65-8005  WEPS  and  U.S.  Air  Force  Contract 
No.  JO  33(615)65-5002.  Results  here  reported  are  preliminary  and  may  bur 
revised  as  additional  data  become  available. 

^Maintained  at  Madison,  Wis.,  in  cooperation  with  the  University  of  Wisconsin* 
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-Underlined  numbers  in  parentheses  refer  to  Literature  Cited  at  the  end  of 
this  report. 
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Much  investigative  work  has  been  done  on  isotropic  cylindrical  shells  sub* 
jected  to  external  pressure  since  Forest  Products  Laboratory  Report  1S44-B, 
giving  theoretical  analysis  of  sandwich  cylinders,  was  published.  It  was 
found  that  experiment  sometimes  yields  critical  loads  that  are  less  than 
those  predicted  by  von  Mises1  theory  (13) .  This  has  been  attributed  to  two 
causes.  First,  the  experimental  cylinders  contained  imperfections  that 
lowered  the  critical  load  (1,  2,  ^3,  <8,  12) .  Second,  energy  levels  associated 
with  postbuckling  configurations  of  the  cylinder  are  lower  than  those  just 
at  buckling.  The  energy  levels  associated  with  postbuckling  may  be  reached 
without  snap- through  buckling.  The  energy  necessary  for  snap- through  buckling 
may  be  supplied  by  vibration  or  shocks  (4,  .5 ,  8) •  The  curves  in  this 

report  do  not  consider  snap- through  buckling  or  cylinders  with  imperfections . 
Sandwich  cylinder  walls,  however,  are  much  more  perfect  than  their  solid 
counterparts  because  they  are  thicker  and  the  effect  of  an  imperfection  is 
in  proportion  to  the  ratio  of  its  amplitude  to  the  thickness  of  the  cylin¬ 
drical  shell  wall.  Also,  the  curves  neglect  the  stiffnesses  of  the  indi¬ 
vidual  facings.  These  stiffnesses  add  to  the  critical  loads  when  the  cylinders 
are  short,  and  it  is  for  short  cylinders  that  snap- through  is  likely  to  occur 
(6). 


Development  of  Formula  for  Design  Curves 


The  formulas  presented  in  this  report  for  calculating  the  critical  external 
pressure  and  the  buckling  coefficient  are  based  on  the  work  presented  in 
Forest  Products  Laboratory  Reports  1844-A  and  1844-B  (9)  with  the  following 
notation: 


Ej.,  E2  Modulus  of  elasticity  of  the  outer  facing  (1)  and  inner  facing  (2), 
respectively. 

Poisson’s  ratio. 


Grz  Modulus  of  rigidity  of  core  in  the  radial  and  axial  directions. 

Gr0  Modulus  of  rigidity  of  core  in  the  radial  and  tangential  directions, 

d  Thickness  of  the  sandwich, 

h  Distance  between  facing  centroids. 


t^,  Thickness  of  outer  and  inner  faces,  respectively, 

r  Mean  radius  of  the  sandwich  cylinder. 


L  '  Length  of  the  cylinder. 
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r^ 


Number  of  half-waves  in  the  circumferences  of  the  cylinder. 
2EjtjE2t2 

3rGr0(l  -  V2)  (E^i  +  E2t2)  * 


Eltl  +  E2c2 


q  The  external  critical  pressure  on  the  cylinder. 

A  few  parameters  of  Report  1844-B  were  generalized  to  dissimilar  facings 
as  follows: 


Report  1844 -B 


Revised  Parameter 


1  +  £  - 


a  Era 


bE^i  E1t-i  leg  b/a 


_  Eea(l  - 


w.  -  E<=a(1  -  »  > 


i  »  ls  2 


a 


L.  i-K?) 

Et 


£2 


These  generalized  parameters  were  combined  with  the  notation  and  derivation 
of  equations  given  in  Report  1844-B.  The  critical  external  pressure  was 
then  found  by  solving  the  revised  form  of  equation  51  of  Report  L844-B  for 
the  buckling  coefficient.  The  determinant  was  simplified  by  taking  modulus 
of  elasticity  of  the  core  (Ec)  to  be  infinite.  For  most  core  materials, 
except  possibly  for  low-density  foams,  Ec  is  sufficiently  large  so  that  this 
assumption  yields  only  slightly  high  values  of  tne  critical  pressure.  Before 
Ec  is  allowed  to  approach  infinity,  the  first,  third,  and  fourth  columns  of 
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the  determinant  are  multiplied  by  SsQ.  Then  when  Ec  approaches  infinity,  the 

Ec  — 

expressions  in  rows  3,  4,  5,  and  6  of  column  3  approach  zero.  Further 
simplification  of  the  6  by  6  determinant  was  made  by  taking  £  to  be  1/3, 
multiplying  row  1  by 


Gz0a(l  ~  t*2) 

E2t2 

and  row  2  by 

Gzea(3-  -  M-2) 

Elcl  b2 

Further  reduction  is  accomplished  by  adding  multiples  of  one  row  to  another 
as  indicated  by  the  following  sequence  of  substitutions,  where  subscripts 
indicate  row  number  counting  down  from  the  top: 


(Rl  +  R2> 


_ EltlE2t2 _ 

Gr0a(l  ~  l1^)  (E^i  +  ^2^2) 


Rj  +  2R2  ■■  R3 

r4  ~  %  - ^  ®4 

%  +  (n2  +  3\2)  R2 - ^  R5 


These  substitutions  cause  the  expressions  in  rows  2,  3,  4,  5,  and  6  of 
column  3  and  those  in  rows  3,  4,  5,  and  6  of  column  6  to  become  zero,  and 
the  6  by  6  determinant  was  readily  reduced  by  minors  to  a  4  by  4  determinant. 
The  radii  a  and  t>  were  eliminated  by  the  following  equations  obtained  from 
the  geometry  of  the  cylinder: 


a  =  r  +  •- 
2 
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and  the  following  substitutions  were  made: 


§  -4s 


After  setting  this  determinant  equal  to  aero  and  solving  for  k,  the  critical 
pressure  can  be  found  by 


Eltl  +  e2c2  , 

q  - 5 —  k 

r(l  -  V-2) 


This  represents  a  theoretical  solution  for  the  critical  pressure  with  the 
assumption  that  the  sandwich  core  modulus  of  elasticity  is  infinite. 

The  determinant  was  further  simplified  by  t  few  approximations  without  any 
significant  loss  of  accuracy.  Examination  of  the  parameters  §  and  4* 

showed  that  §  is  large  in  comparison  to  4?.  Therefore,  any  terms  of 
(§  +  4*)  were  taken  as  equal  to  $  . 

Making  this  simplification,  the  determinant  reduces  to  the  following  ex¬ 
pressions: 

Sow  1,  column  1 

- - 2 - +  2  y  (l - 

(R+l)n2  2  V  \  3ny 

Bow  2,  column  1 

-4+  6V 

n* 


4>. 

T 


Bow  3,  column  1 

— 1?-  <n  +  3»4-  +  3vp 

(R+l)2  n2  i 


Bow  4,  column  1 


8R  (n2  -  D 
(R  +  l)2  n2 


4> 

(n2  +  3p)-|-  +  3Vp 
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Row  1,  column  2 


+  |+  <n2  -  1)  (2k  -  1) 


Bov  2,  column  2 


Row  3,  column  2 


20 


sL^-k  +  fj  +  <n2  -  1)1^  +  2k(n2  -  lHn2  +  30) 


Row  4,  column  2 
/ 


20 


^  fcLz.il  +  p  +  (n2  -  1)1^  +  2k(n2  -  l)(n2  + 


30) 


Row  1,  column  3 


2  4> 

§ 


Row  2,  column  3 
0 


Row  3,  column  3 


.  *±JL  +  („2  +  3p)^fc 


Row  4,  column  3 


R  +  1  +  (n2  +  30)M 


Row  1,  column  4 


n2  £ 

3 
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Row  2,  column  4 


4 

3  § 


Row  3,  column  4 


*3n2  ^  -  (n2  +  3p)|  0 

Row  4,  column  4 

40n2  ^  -  |V  <n2  +  3£)  j  p 


By  setting  this  determinant  equal  to  zero,  the  approximate  solution  for  k  is 

.. .  I  *  Trfrr?  -»{3+ r)[(r  ~  l)<°2  - 1  +  » -  I  42  + 1] 

[(r + 1)2<“2  • 1}  +  ^ +  3w,(n2  ^  f)] 

where  a  *  JL.  *  ~ 

$  2r 

2  2 

To  obtain  this  value  of  k,  terms  containing  k__  and  kq  were  neglected, 
was  also  assumed  that  terms  (1  ±  ma)  =  1  where  m  is  a  small  whole  number 

A  loner  and  upper  limit  exist  for  the  buckling  coefficient  k.  The  lower 
limit  is  associated  with  the  infinitely  long  shell  for  which  p  =  0  and 
this  limit  from  the  approximate  formula  for  k  is  given  by 

p  ,  0;  k  - - £Rfc?  -  1)q2 

(R  +  l)z(l  +  3VOnz) 


and  for  V*  0  this  formula  is  minimum  for  n  =  2  and  is 

12  S32 


P  a  0;  V  *  Oj 


(R  +  l)2 
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The  upper  limit  for  ^  is  associated  with  the  usual  shear  instability  type  of 
buckling  of  the  sandwich  wall  which  occurs  when  n  =  «  and  the  formula  for  k 
then  is  given  by 


n  a  ea  ; 


k’BCL 

3(R  +  l)2v 


Substitution  of  this  value  for  k  into  the  formula  for  q  results  in  a  critical 
hoop  compression  per  unit  length  of  cylinder  of  N  =  qr  =  hG . 


Description  of  Design  Curves 


Using  the  approximate  equation,  values  for  k  were  determined  for  various 
values  of  V,  n,  i,  and  a£.  These  values  of  k  are  plotted  (figs.  1  to  16)  for 
values  of  —  ranging  from  1  to  100.  Curves  are  gi/en  for  values  of  V  equal 


to  0,  0.5,  1.0,  and  1.5.  For  each  of  these  valu  .s  of  V,  curves  are  given  for 

seven  values  of  o£  ranging  from  10"  ^  to  10*^.  curves  apply  to  sandwich 

cylinders  with  isotropic  facings  and  isotropic  cores  or  orthotropic  cores 
with  their  natural  a;:es  parallel  to  the  axial,  tangential,  and  radial  direc¬ 
tions  of  the  cylinders.  The  critical  pressure  is  given  by 


q 


k 


From  a  solution  of  the  exact  determinant ,  it  was  found  that  the  modulus  of 
rigidity  of  the  core  associated  with  the  radial  and  axial  directions  Gj^ 

has  very  little  influence  on  the  critical  pressure.  It  does  not  enter  the 

GtJ_ 

formulas  for  the  limits  of  critical  pressure.  The  ratio  6  =  -JSs  does 
affect  the  critical  pressure  for  small  values,  but  for  reasonable  values 


of  jB  the  critical  pressure  is  affected  onl^  slightly.  The  ratio  B  does  not 
appear  in  the  approximate  solution  for  k. 


The  accuracy  of  the  approximate  solution  was  checked  by  making  a  comparison 
of  the  approximate  and  exact  solutions  for  values  of  V  =  0  and  (X  -  0.01. 
This  comparison  is  shown  on  the  top  curve  of  figure  1  where  the  small 
circles  represent  the  exact  solution.  This  shows  that,  except  for  small  i 

values,  the  error  encountered  with  the  approximate  solution  is  less  than 
1  percent .  For  larger  values  of  V  and  smaller  values  of  Of,  the  error 

diminishes.  For  small  values  of—,  it  is  doubtful  that  either  solution 

r 

gives  good  results.  — 
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Figure  13. --Values  of  k  for  V  *  1.5,  and  for 
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Figure  16,— Values  of  k  for  V  »  1.5,  and  for 
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pressure  of  finite  length,  circular  cylindrical  shells  pressure  of  finite  length,  circular  cylindrical  shells 

with  sandwich  walls.  with  sandwich  walls. 


FOREST  PRODUCTS  LABORATORY 


The  following  lists  of  publications  deal  with  investigative  projects  of  the 
Forest  Products  Laboratory  or  relate  to  special  interest  groups  and  are  avail¬ 
able  upon  request: 


Architects,  Builders,  Engineers, 
and  Retail  Lumbermen 

Box,  Crate,  and  Packaging  Data 

Chemistry  of  Wood 

Drying  of  Wood 

Fire  Protection 

Fungus  and  Insect  Defects  in 
Forest  Products 

Furniture  Manufacturers, 
Woodworkers,  and  Teachers 
of  Woodshop  Practice 

Glue  and  Plywood 


Growth,  Structure,  and 
Identification  of  Wood 

Logging,  Milling,  and  Utilization 
of  Timber  Products 

Mechanical  Properties  of  Timber 

Structural  Sandwich,  Plastic 
Laminates,  and  Wood-Base 
Components 

Thermal  Properties  of  Wood 
Wood  Fiber  Products 
Wood  Finishing  Subjects 
Wood  Preservation 


Note:  Since  Forest  Products  Laboratory  publications  are  so  varied  in  subject 
matter,  no  single  catalog  of  titles  is  issued.  Instead,  a  listing  is  made  for 
each  area  of  Laboratory  research.  Twice  a  year,  January  1  and  July  i, 
a  list  is  compiled  showing  new  reports  for  the  previous  6  months. 
This  is  the  only  item  sent  regularly  to  the  Laboratory’s  mailing  roster, 
and  it  serves  to  keep  current  the  various  subject  matter  listings.  Names 
may  be  added  to  the  mailing  roster  upon  request. 


